Cause of Disparity between the Results for the nn 

Transition Problem 

V.I.Nazaruk 

Institute for Nuclear Research of RAS, 60th October 
Anniversary Prospect 7a, 117312 Moscow, Russia 
E-mail: nazaruk@al20.inr.troitsk.ru 

Abstract 

We show that there is a double counting in the standard model of nn mixing in the medium, 
resulting in full cancellation of leading terms. The direct calculation of nn transition, annihi- 
lation is performed. For lower limit on the free-space nn oscillation time we get r^jn ~ lO^^y. 
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Any information on the occurrence of nn oscillation [1,2] is important in order to discriminate 
among various grand unified theories. The most direct limit on the free-space nn oscillation time 
Tnn is obtained using free neutrons: Tnn > 10^s[3]. Alternatively, a limit can be extracted from the 
nuclear annihilation lifetime T measured in proton-decay type experiments: Tnn > lO^s ~ ly (see, 
for example, Ref. [4]). The calculations involved were based on the potential model of nn mixing in 
the medium. In this letter the model independent approach is presented. 

l.In Ref. [5] for free-space nn oscillation time the limit T^in — 3 • lO^^y was obtained, which 
increases the previous one (see, for example, Ref. [4]) by 31 orders of magnitude. First off all we 
expose a drawback hidden in the standard model. In the standard approach (labelled bellow as 
potential model) the nn transitions in a medium are described by Schrodinger equations 

{idt + VV2m - Un)n{x) = en{x), {idt + W'^ /2m - Un)n{x) = en{x). (1) 

Here e = I /Tnn is a small parameter [4]; Un and Un are the self-consistent neutron potential and 
n-nucleus optical potential respectively. For Un — const, and Un — const, in the lowest order on e 
the probability of the process is 

W.,ot{t) = 1- I Uu{t) P= 2ImTu{t), Tu{t) = t{e/6U)^[l - t6Ut - exp{-t6Ut)], (2) 

where U — 1 + iT, T are the evolution and T-operator respectively; 

5U^Un-Un, Un^ReUn-ir/2, (3) 

r ~ lOOMeV is the annihilation width of n-nucleus state. 

What is meant by Wpot{t)7 Let us take the imaginary part of Eq.(16) of Ref.[5] 

2ImTii{t) = eH^ - /* dta t" dt82ImT^i{T), (4) 

Jo Jo 

T = ta—tp. Here T{t) , T"'{t) are the T-matrix of the whole process and n-nucleus decay respectively. 
From the condition of probability conservation Z^/ | 1 we have 

2ImTu = E I I' + I I'' E I ^/iW l'= ^W- (5) 

m m 

The process probability W{t) will be specified bellow. | T/j p~ e^, whereas | 1^6. So for 
the l.h.s. of Eq.(4) 2ImTu{t) = W{t), that was taken into account in (2). For the T-matrix of 
annihilation nucleus decay J^"(t) Eq.(5) has the form 

2/m7^(r) = E I T%{t) r + I nir) ^ . (6) 

The annihilation nucleus decay is nonperturbative process and | S/^i I ^/i Now Eq.(4) 

has the form 

W{t) = eH' - f dta f" dtp I T^,{t) \^ -e^ f dt^ dtp E I Tf,{T) \' (7) 



Let us calculate and Tj* in the framework of potential model. The wave function of initial state 
obeys equation 

i^^Ho<^, Ho = -Vy2m + Un. (8) 
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In t = the interaction SU is turned on. We have 



.9* 



{Ho + 5U)'i/, ^(0) = $(0). 



(9) 




The projection to the initial state and T-matrix at t — t are 



< $ 



*>=C/^(r)=expH5[/r), 



(10) 



T^{t) = t[l - exp(-z5t/r)], ^ | T^ir) \'= 1- | U^{t) \'= 1 - e"^- = W^,(r), (11) 



where Wn{T) is the n-nucleus decay probabihty. Note that F corresponds to all n-nucleus interactions 
followed by annihilation. However, the main contribution gives the annihilation without rescattering 
of n[6], because aann > '^cTsc- Substituting these expressions in (7), one obtains the potential model 
result (2). 

Therefore, the finite time approach was verified by the example of exactly solvable potential 
model. It is involved in Eq. (4) as a special case. 

Let us return to Eq.(7). It is at least unclear. l.The first term is free-space nn transition 
probability. Matrix elements T-^ and T"_^,- describe transitions {n — nucleus) (n — nucleus) and 
(n — nucleus) {annihilation products) respectively. So the first and second terms correspond to 
n-nucleus in the final states. However, in the experiment only annihilation products are detected and 
the result should be expressed as T^^j solely. Moreover, n-nucleus decays into final state products 
identical with those given by third term. This suggests that potential model contains the double 
counting. Expression 1— | Uu p from Eq.(2) describes the inclusive decay of initial state and so 
the nn transition with n-nucleus in the final state is also included in Wpot, unless additional limits 
are imposed. To obviate the double counting the annihilation products in the final state should be 
fixed. 2. Let us | 6Ut \<^ 1. (This regime occurs in other problems.) When F = 0, the third term 
equals to zero. When F 7^ 0, the contribution of the third term is negative and dW/dT < 0, whereas 
the opening of the new channel (annihilation) should increase W. 

How much is the probable error? Contributions of the second and third terms arc: X2 = 
—e^t"^ /2 + F2, xs = —e^t^/2 + F^. Functions ^2,3 contain the terms proportional to t and exp(— i5t/i). 
So the eH"^ term produced by third term is fully canceled. This is a consequence of double counting. 
Therein hes a reason of discrepancy between ours and potential model results. Solving Eqs.(l) by 
method of Green functions we will reach the same results. We have started from Eq.(4) only for 
verification of finite time approach. 

As noted in[5], Eqs.(ll) and (2) can be also obtained by means of microscopic variant of potential 
model (zero angle rescattering diagrams of n) . In this case the Hamiltonian of n- medium interaction 
is H — 5U. The same calculation was repeated by Dover et al.[4]. They substitute H — — iF/2 in 
(4) and obtaine (2). On the basis of this and only this they refute the result of Ref. [5]. In other 
words they refute our limit because it differs from prediction of potential model {H = —iF/2 [4]). 

Our concern is with Z^/^j | Tfi p. It is connected with diagonal matrix element by Eq.(5): 



Calculation of Tu is determined by r.h.s., namely, the cut corresponding to Ta must contain only 
annihilation products, that is not in accordance with Eq.(7). It includes superfiuous "incorrect" 
states /' = {n — nucleus) . Note, that eigenfunctions of Hq + 5U do not form the complete 



(12) 
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orthogonal set. Because of this the {n— nucleus) (described by Un) also may not appear in Eq.(12) as 
intermediate state. So model (1) is inapplicable to the problem under study because it automatically 
leads to incorrect matrix element Ta. Elimination of superfluous trajectories from Ta means the 
direct calculation of T/,. We can write W — 2ImTii, however Ta should be calculated by means of 
r.h.s. ofEq.(12). 

2. In Rcf. [5] the first and third terms were taken into account. The second one was omitted. The 
first term reproduces low density limit and has a sense for nn transitions in the gas. This scheme 
is not quite correct. In this paper we present the direct calculation of the process {nucleus) 
{n — nucleus) — > {annihilation products) . We have 



< / I f/(t,0) - / I Onp >=iTfi{t) = 
Er=i(-^)'+' < / I /oC^ti-./o^'c^^fc/o' dtpH{ti)...H{tk)Hr,n{tp) I Onp >, 



(13) 



where 

H{t) — {all n — medium interactions) — Un, Hnn{t) — e j d^x{^n^n + h.c), (14) 

H + Hnn = Hi. Here | On^ > is the state of medium containing the neutron with 4-momenta 
P = + Un, Pn), < / | IS the annihilation products; Hnn is the oscillation Hamiltonian[4]. In 

the case of the formulation of the (S-matrix problem (t, 0) (oo, — oo) Eq.(13) in the momentum 
representation includes the singular propagator G = l/(e„ — p^j/2m — Un) ~ 1/0. Taking into 
account that Hnn |Onp>= e |Onp>, we change the integrating order and obtain 

Tf,{t) = -eJ^dtpiTf,{t-tp), 
iT%{T) = Er=i(-^)'/i dt^.-.f,;-' dtk < f I H{t,)...H{t,) I On, >, ^'^> 

where T^^ is an exact amplitudes of n-nucleus decay, | Oup > is the state of medium containing the 
n with 4-momenta p; t — t — tp. 4-momenta of n and n are equal. 

The 2-step process was reduced to the annihilation decay of n-nuclcus. (The slightly different 
method is that antincutron Green function is separated [5].) It is seen from (13), and (15) that 
both pre- and post- nn conversion spatial wave function of the system coincide 

\Onp>sp^\Onp>sp . (16) 

We would like to stress that in potential model the picture of n- nucleus formation is precisely the 
same. Really, let us the nn conversion takes place at t = 0. Solution of Eqs.(l) is continuous and 
\[/(-0) = n = \&(+0) = n, that is identical to (16). (Sec also Eqs.(8),(9).) Note that Eq.(4) was 
obtained in perfect analogy to (15). In particular, for and condition (16) was fulfilled. Here- 
after, the potential model of n-medium interaction (block T") was used and Wpot was reproduced, 
which also corroborates the picture of n- nucleus formation given above. 

In both models the stage of nn conversion is identical. The basic difference centers on the next 
stage - annihilation. In the potential model is calculated end used in Eq.(7), which is wrong. 
We calculate T^j starting from the same point (16). The result will be expressed through T (see 
Eqs.(19),(ll)), but not 5U, as it usually is in calculation of decays. 

Let us take t — T — ■ 10^^y[7], where T is nuclear annihilation hfetime. The characteristic 
annihilation time of n in nuclear matter is l/F ~ 10~^^s. When r » l/F, Tf-{T) reaches its 
asimptotic value T^f 

T%{t » 1/r) = T;(oo) = T% = const. (17) 
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The expressions of this type are the basis for all S'-matrix calculations. (Measurement of any process 
corresponds to some interval r. Consequently, it is necessary to calculate U{t). Replacement 
U{t) — > S{oo) is equivalent to (17).) From (15) and (17) we have 



Tf,{T) = -^e[ | dt^Tf.iT - tp) + j dtpT%{T - tp)] ~ -ieTT%. (18) 



rT-l/T _^ _ rT 

h-i/r 

Really, | Tf^(r) |< 1 because l^/i(r) =| T^t) p< 1; T = 2 ■ lO^^^s, or T ~ lO^Vr in units of l/F. 
Obviously, the contribution of second term is negligible. The probability of the whole process is 



W{T) = Yl I ^/i(^) ^^^^ Yl I ^/i(^) ^^'^^^ (19) 

where Eq.(ll) have been taken into account. The limit for is obtained from the inequality 
W{T) < 1. For T = 6.5 • 10=^^y[7] we have Tmin ~ lO^^y. 

3. Let us return to the reason of enormous quantitative disagreement between the our and poten- 
tial model results. The strong result sensitivity was to be expected. Really, the ^-matrix amplitude 
Ms, corresponding to nn transition, annihilation diverges: 

Ms = e ^ —M ~ 1, (20) 

where M is the annihilation amplitude. This is infrared singularities conditioned by zero momentum 
transfer in the e- vertex. It is easy to understand that Mg ~ 1/0 for any bound state wave function of 
neutron (i.e., for any nuclear model). On the other hand from Eqs.(l) it is clear that in the potential 
model the energy is not conserved and becomes complex in the e- vertex Ma — > Ma + SU {Ma is the 
nuclear mass). The corresponding antineutron Green function is G = l/(e„ — p^/2m — C4) = 1/5U. 
SU = is the peculiar point of Mg. So Mg is extremely sensitive to 6U. (Usually, the SU- dependence 
of G is masked by momentum transferred q: = (e„ — qo) — (p„ — qf' /2m — Un — SU. We deal 
with 2-tail and q = 0.) 

Comparing (20) with (18) one sees that primitively the limit 5U — > corresponds to replacement 

1/SU^T. (21) 

Certainly, we do not set 5U — 0, because S'-matrix amplitude is not considered at all. In the 
calculation of Eq.(13) the multiplier T (see Eq. (18)) arises automatically instead of 1/Ag in the 
S'-matrix theory. When g ^ in the e- vertex, Eq.(13) leads to usual S-matrix result (see below). 
Formal reason for the differences in the results is the full cancellation of terms ~ in Eq.(7). 
Erroneous structure of (7) conditioned by nonperturbative and 2-step character of the process. 
q — extremely reinforces the disagreement. 

4. One additional comment is necessary regarding t-dcpendence of the whole process probability 
W{t). Eq.(19) has been obtained in the lowest order on e. The precise distribution Wprit) which 
allows for the all orders on e is unknown. However, W is the first term of the expansion of Wpr and 
we can restrict ourselves to a lowest order Wpr = W, as it usually is for rare decays. Wpot is also 
calculated in the lowest order on e. 

The protons must be in very early stage of the decay process. Thus the realistic possibility 
is considercd[8-ll] that the proton has not yet entered the exponential stage of its decay but is, 
instead, subject to non-exponential behavior which is rigorously demanded by quantum theory for 
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sufficiently early times. At ffist sight, since Tnn > lO'^^y for nn-mixing in nuclear the similar picture 
should be expected. In fact the situation is more serious. We deal with two-step process. In 
attempting to calculate Mg and ~| ^ in the framework of standard S- matrix theory we get 

~ 1/0. So there is no sense to speak about decay law exp(— Fst). It is necessary to calculate the 
distribution W{t) as it was done above. 

Finally, we will touch upon the main points of Krivoruchenko's preprint [12]. (1) The nn tran- 
sition, annihilation (two-step nuclear decay) and particles motion in the classical fields are the 
different problems. Describing the first one by Eqs.(l) we understand that this is an effective pro- 
cedure. Prom formal standpoint in the first and second cases the potentials are complex and real 
respectively. Unfortunately, sometimes the literal analogy between these problems is drawn[12]. (2) 
The initial Eq.(ll) of Ref. [12] must describe the nn transition, annihilation. However, the l.h. s. 
of Eq.(ll) is free of n -nucleus interaction at all. The r.h.s. contains annihilation width F (we stress 
this point) and coincides with potential model result. We would hke also to get the result without 
calculation, but some difficulties emerge in reaching this goal. 

We attempt to calculate the process amplitude starting from (14). The S- matrix theory gives 
(20). The approach with finite time interval is infrared- free. Its verification for diagrams with q = 
was made above by the example of potential model. For nonsingular diagrams the test is obvious. 
Let us g 7^ in the e-vertex. Appropriate calculation with finite time interval (adiabatic hypothesis 
should be used) converts to the S'-matrix result: T/, = ie'{l/ Aq)Tjl, where Tjl is the -S-matrix 
amplitude of annihilation of virtual n with 4-momenta k = p — q. Comparing with (18) one sees 
that limit Aq corresponds to replacement 1/Ag t (compare with (21)). Similar problem 
for matrix element was solved in Ref. [13]. Note, however, that there are essential differences 
between above mentioned problems. This result as well as connection between jS-matrix theory and 
approach with finite time interval for the diagrams of various type will be presented in the next 
paper. We will also considered the another exactly solvable problems and show that all the results 
are true for any nuclear model. 

In conclusion, we perform the direct calculation of the process amplitude with annihilation 
products in the final state. Potential model is inapplicable to the problem under study. This 
explains the different functional structure of the results: W{T)/Wpat{T) ~ FT. If it is remembered 
that T ~ lO^^y the quantitative distinction becomes clear as well. 
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